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ABSTRACT 

Estimating the number of microlensing events observed in different parts of 
the Galactic bulge is a crucial point in planning microlensing experiments. 
Reliable estimates are especially important if observing resources are scarce, 
as is the case for space missions: K2, WFIRST, and Euclid. Here we show that 
the number of detected events can be reliably estimated based on statistics of 
stars observed in targeted fields. The statistics can be estimated relatively eas¬ 
ily, which makes presented method suitable for planning future microlensing 
experiments. 
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1 INTRODUCTION 


The scientific yield of any gravitational microlensing experiment strongly depends on the 
number of detected events. A higher number of events improves the statistical properties 
of the event sample and allows one to detect larger number of microlensing anomalies, e.g., 
planetary signals provided that the observing cadence is high enough. Hence, it is of primary 
importance to predict how many events would be detected in a given part of the Galactic 
bulge . Additionally, the event rate gives us information about the structure of the bulge 


(e.g., 


Kiraga fc Paczyhskilll994l ). The connection betw een statistics of observed events and 


the structure of the bulge is an area of active research flKerins et ah 


2009 


Snmi et ahl 


2013 
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Wvrzvkowski et al. 


2015h . but a full understanding of microlensing event rate has not yet 


been achieved. 

Here we present a new approach to a problem of predicting microlensing event rates. 
We show that the number of events observed scales with the product of the surface density 
of stars brighter than the completeness limit (W(/ < Rm)) and the surface density of 
red clump (RC) stars in given held (A^rc)- The RC stars are relatively bright and easy to 
identify. Hence, the RC sample is essentially complete in any given held, except the highest 
extinction helds. Thanks to their brightness, number of RC stars is a good proxy for total 
number of stars observed in given sight-line i.e., they estimate the number of potential lenses 
in given held. The surface density N^,{I < RR) is the proxy for a number of potential sources 
in microlensing events. Hence, one may conjecture that the number of events observed is 
proportional to the product of number of potential lenses and number of potential sources. 
This is obviously true if the stellar kinematics and the distribution of stars along the line 
of sight did not changed from held to held. It is important to note that both < RjR 
and A^rc can be relatively easily measured. The next section summarizes the data used. In 
Section 3 we present the basic correlation we found. Section 4 describes ehorts to hnd tighter 
relation. We discuss results in Section 5. 


2 DATA 


We an alyse the catalogue of standard microlensing events presented by 


Wvrzvkowski et al. 


(1201511 . The events are standard in the sense that they do not show obvious signatures of 
either binary lenses or sources, nor extended source ehects. Microlensing parallax ehect may 
be pronounced in some of these events but they are considered standard as long as the light 
curve is well htted by a point-source point-lens model. We prefer a catalogue of standard 
events selected after the observations are completed over the list of candidate events selected 
while the observations are ongoing, because the former is more objective, more unifo rm, more 


complete, and contains smaller number of false-positives. 


Wvrzvkowski et al. 


(120151) searched 


photometric databases of the third phase of t he Optical Gravitational Lensing Experiment 


(OGLE-HI; 


Udalski 

20Q3; 

Udalski et al. 

2ood 


20081 ). The sky images were collected between June 


2001 and May 2009 using the 1.3-m telescope at Las Gampanas Observatory, Ghile and an 
eight chip camera with a held of view of 36' x 36'. Microlensing events were searched after 
the project was completed, thus events peaking during whole project time span were found. 
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We note that 15 percent of sources are RC stars. Different astrophysical phenomena that 
can mimic microlensing events, like dwarf nova outbursts or brightening of Be stars, were 
rejected thanks to the long term photometry. 

There are a variety of factors that affect the number of detected microlensing events, 
some of which are of astrophysical origin and other are survey-specific biases. An important 
factor, that we are not able to fully eliminate, is different observing cadence in different 
fields. Additionally, the number of epochs in a given field changes from one year to another. 
Differing observing cadence affects the efficiency with which surveys can hnd events, espe¬ 
cially short events. We minimize the effect of changing cadence by selecting for our analysis 
a subset of fields and a subset of observing seasons. We remove events peaking in 2001, 
as the OGLE-III survey started in the middle of the 2001 bulge season and the observing 
strategy differed significantly from other years. The last year of the OGLE-III observations 
ended before the best observing conditions for bulge helds started in June. Thus, we also 
eliminate events peaking in 2009. Other than that, the biggest change in observing strat¬ 
egy was introduced in 2004 when many helds were abandoned and the observations focused 
on helds with higher number of events; hence, we decided to analyse only events peaking 
between 2004 and 2008. Even after imposing these criteria we are left with helds observed 
at diherent cadence. In order to impose reasonable uniformity we excluded helds with less 
than 90 epochs per year. This leaves us with 63 helds for which mean number of epochs per 
year varies between 94 and 262 with average of 165 and a total of 2563 events. The centres 
of these helds are within Galactic coordinates range —4?9 < I < 4?6 and —4?8 < b < —1?3 
and are outlined in Fig. [H 

In addition to the information about events, we also need data that describe the structure 
of the Galaxy. We employ two kinds of data: the surface density of RG stars and the surface 


density of all stars down to a certain 


investigated in detail by 


Nataf et af 


i mit. T he RG stars in the OGLE- III bulge helds were 


(120131 ) and we use their results. 


Nataf et al. 


(120131) 


divided the OGLE-III helds into a number of sight lines and provide spatial density of RG 


stars for each sight line sepa rately. In the area considered here, three helds have subheld: 


for which 


Nataf et al. 


(j201^ could not measure the properties of RG because of extremely 




high extinction. These sky-areas are not used in our study. We divided every subheld into 32 
regions and interpolated RG spatial density values to the centres of the regions. To hnd the 


^ Subfield is sky area corresponding to single CCD chip i.e., 1/8 of camera field of view. 
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Figure 1. Analysed OGLE-III fields (color-shaded squares) and microlensing events (black points) in Galactic coordinates. 
Fields are divided into three groups marked by different colors as described in Sec. [5] Gray shaded fields are not taken into 
account because of too small number of epochs. Three red-shaded fields have subfields with extinction that is too high to 
measure RG stars and hence we omit them in present study. 


A^rc for e ach subfield we ave raged the interpolated values. In this way we cope with varying 

(120131) sightlines and edge effects. 


density of 


Nataf et ah 


The last piece of information we need is spatial density of stars o 


area. We use the photometric catalogues presented by 


Szvmahski et al. 


□serve d in each sky 


( 2011 1 that list stars 


detected on the OGLE-III reference images. In each subheld we count the number of stars 
down to spe cihed brightness l imit a nd scale it to 1.0 deg^. Note, that each reference image 


analysed by 


Szvmahski et al. 


( 120111 ) was constructed using a set of high-quality images, but 


slightly shifted with respect to a dehned centre. This renders the outer parts of the subhelds 
of lower photometric quality and hence they are not considered for estimating spatial density 
of stars. 


3 PRIMARY CORRELATION 


For each held we derive mean number of events per year per deg^ ( 7 ; this symbol is used 
in analogy to frequently used T as an event rate corrected for the detection efficiency). We 
dehne microlensing event sensitivity as: 


5 , 




a,him — 


RC 


10^ deg“^ 


X 


A,(/< Ju^) ' 
10® deg“2 


( 1 ) 
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We start with a simple choice of a = 1 and Jum = 20 mag. This limiting brightness is 
chosen because it is a rough estimate of the OGLE-III completeness in the dense stellar 
helds. For reference we provide the average numbers for the well-studied Baade window 
(/ = 1?0, b = —3?9): iV^c = 49.9, N^{I < 20mag) = 4.05, and < 20.5mag) = 4.68 in 
the same units as in Eq. [H We present plot of 7 as a function of S'i_ 2 o iii Fig- El The line 
htted to these data is: 

--1-r = 0.0999(44) X ^ 1,20 - 3.21(97) (2) 

deg ^yr ^ 

with x^/dof = 1.46. The uncertainties in S' 1^20 are negligible compared to the Poisson noise 
of 7. It turns out that even our simple approach for Ending sensitivity results in a relatively 
tight correlation. 


4 DETAIL PARAMETER SELECTION 


There are multiple efiects that may afiect the number of potential lenses e.g., the number 
of disk to bulge stars changes with Galactic latitude. Yet, our definition of sensitivity takes 
only Arc into account. Hence we allowed a to vary and selected the value that gives the 
tightest relation for given 

We further tried to improve our event rate estimates by revising microlesning events 
included. In the fields with smaller number of epochs the shortest events are strongly affected 
by selection biases, thus, we checked the impact of rejecting all events with Einstein time- 
scale shorter than some limit tE,iim (set to between 5 d and 10 d). The reduction of the scatter 
in the 7(fE > relation is only slight. In our final calculation we used tE,iim = 8 d, 

which resulted in removing 11 per cent of the event sample. In eighteen best-observed fields 
(at least 1200 epochs in OGLE-III) 14 per cent of the sample was removed. 

An important simplification we made in previous section is that our proxy for number 
of potential sources is the number of stars down to a sharp brightness limit Jiim = 20 mag. 
In a real experiment, no such sharp limit exists; sources can be significantly fainter than the 
photometric detection limit and still be detected if they are highly magnified. On the other 
hand, almost all events with bright sources are easily found. Instead of detailed modelling 
of detection efficiency that could take into account the event time-scales, observing cadence, 
luminosity function etc., we try to improve our relation while retaining its simplici ty. We 


calculate the source brightness (J^) for events presented by 


Wvrzvkowski et al. 


( 2015 ) based 


on their blending ratios and baseline brightness. In Fig [3] we show a histogram of J*. The 
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Figure 2. Relation between number of observed microlensing events in a given field and product of surface densities of stars 
brighter than 20 mag and RC stars. The single point with large errorbar represents field blgl94 in which only a single subfield 
has A^rc derived. 


number of detected events falls to half of maximum at / = 20.5 mag, and this limit is used 
for further calculations of number of stars. 

Our dual relation is presented in Fig. HJand fitted with; 

The resulting x^/dof = 1.17 and is significantly smaller than in primary relation. Fig. IHalso 
presents how residuals from htted relation depend on Galactic longitude and latitude. No 
significant dependence is seen. This shows robustness of our method. 
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Figure 3. Histogram of source 


brightness in IWvrzvkowski et alJ ll2015l^ sample. The bins are 0.2 mag wide. 


The two relations presented in Eq. [2] and [3] show that the choice of a and Jum may 
significantly impact how tight the relation is. In order to fnrther verify how optimal valne 
of a depends on Jum we present a plot of y^/dof as a fnnction of Jum and a. One clearly 
sees that both parameters are anti-correlated and the best hts correspond to Inm in a range 
from 19.9 to 20.6 mag. 
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Figure 4. Final relation between number of observed microlensing events with t-E > S d in given field and -So. 55 , 20 .5 sensitivity. 
Right panels present residuals as a function of Galactic coordinates. 


5 DISCUSSION 


We have shown that the number of observed microlensing events can be estimated based on 
the statistics of stars observed in a given held. It is important to note that presented frame¬ 
work takes into account all stellar populations that constitute the bulge and disk, including 
the ones that have not yet been fully described. As an example, the X-shaped structure was 


hidde n for a long time and discovered only recently flNataf et ah 


2010 


McWilliam fc Zoccali 


2 OIOI) . Our method takes such structures into account via N^,{I < I\[^) and, unlike efforts 
that depend on the Galactic distribution models, does not depend on current knowledge of 
different stellar populations. 

We note that better understanding of different factors contributing to bulge microlensing 
could be achieved via a two-step investigation. First, the Galactic distribution models would 
give reliable predictions of the number density of brightness-limited sample of stars stars and 
the number density of RG stars. Second, one would use relations presented here to estimate 
the expected event rates and compare these to the measured values. The events that appear 
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Figure 5. The x^/dof surface in — a space. 


in the sky-areas with low sensitivity have higher chance of being disk-disk events. There 
is particular interest in events involving globular cluster lenses. Cluster lenses h ave known 


dista i ices and proper motions, hence one can measure their mass based on t^ only fjPaczyhski 


1994 


Pietrukowicz et al 


] 20121 . 


Our framework could be used to select promising candidates 


for cluster lensing events: the events appearing in low sensitivity regions and near the central 
regions of a globular cluster. 

The key application of the relations presented above will be in selection of particular 
helds for satellite survey observations. The planning of survey helds requires comparison of 
expected yields on different angular scales. We tried to verify if our relations hold in different 
areas of Galactic bulge in two different ways. First, we divided analysed helds into three 
groups with following constrains: each group represents continuous (or near-continuous) sky- 
area, the three helds with the highest event rate are in diherent groups, and total number 
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Color in Fig[T] 

^fields 

^events 

< 1 > 

<b> 

a [deg 2yr ^ 

b [deg 2yr 1] 

X^/dof 

blue 

18 

853 

3? 07 

-2? 86 

0.792(76) 

-15.3(33) 

2.11 

green 

20 

855 

0?90 

-3? 45 

0.734(56) 

-12.7(23) 

0.62 

red 

25 

855 

-2? 40 

-2? 92 

0.774(65) 

-15.6(26) 

1.00 


Table 1. Comparison of event rate vs. sensitivity relation fits in groups of fields outlined in Fig^by different colours. Columns 
give: colour of shading in Fig^ number of fields in a group, total number of events, mean Galactic coordinates, and linear fits 
results in the last three columns. The relations fitted were in the form: 7(tE > 8d) = a x 50.55,20.5 + be., the same as in 
Eq.El for whole sample. 


of events in each group is similar. The division of OGLE-III helds in the three groups is 
presented by different colours in [H The groups differ mostly in mean Galactic longitude. 
The parameters of each group and the results of line htting are presented in Table [T] There 
are no signihcant differences in values of parameters htted between different groups as well 
as between them and relation presented in Eq. [3l On the other hand, the y^/dof varies by 
more than a factor of three. It is the largest in the East part of bulge. 

The above test verihed how our empirical relation works on scales of a few degrees, but 
choosing optimal helds for microlensing observations requires also a comparison of sky-areas 
on arcminute scales. Thus, the second way we tested our relation is based on 5 arcmin^ sky- 
regions used to estimate A^rc- For each region we calculated sensitivity S'0.55,20.5 and expected 
event rate. In Figure |6] we compare the measured and expected event rate for all regions 
with sensitivity S'0,55,20.5 < *S'o,55,2o.5- Note that to hnd expected event rate we integrate over 
small regions that are well separated on the sky. The measured and predicted values are very 
similar for the range of S'0,55,20.5 used in deriving Eq. [3]i.e., for S'0.55,20.5 > 25 the relative 
difference is smaller than 3%. On the other hand, for smaller So.55,20.5 fhe differences are 
larger e.g., for S'0.55,20.5 = 12.5 the relative difference is 21%. The overall similarity of the 
two curves in Fig. |6] proves that our method reliably predicts event rate. 

The empirical relation given in Eq. [2] and |3] was derived based on a subset of OGLE-III 
bulge helds but can be used to predict the event rates in other helds as well. We estimate the 
values A^rc and < 20.5mag) in the same way as previously and average the expected 
event rate 'y{t-E > 8 d) over the area of a subheld. For 2% of the subhelds the A^rc values 
are not available. The expected event rates are presented in Figure [71 The highest value of 
7 (fE > 8 d) are close to 50deg“^yr“^ and are concentrated in small sky area I ~ 0?5 and 
b ^ —2°. Most of the sky-area at |/| > 6° or |6| > 5° shows very small event rates. All the 
estiamted values are provided in Table |2l 

Our empirical relations are based on a subsample of helds and events observed by the 
OGLE-III survey. Similar relations can be found for other high-cadence bulge observing 
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Figure 6. Comparison of predicted and measured event rates for all 5 arcmin^ sky-regions with sensitivity 50.55,20.5 smaller 
than given value *§ 0 . 55 , 20 .5- The predicted event rate is multiplied by 1.03 to match the total observed event rate. This correction 
is smaller than the uncertainties in Eq. [3] 


surveys if corrected for different cadence, band-pass etc. The main application of such efforts 
should be selection of helds for microlensing surveys. We note that held selection does not 
depend on the zero point of scaling relation; it is enough to hnd relative number of events 


that a given survey should hnd in diherent helds to select optima 


highest pri ority i.e., space-base d surveys, are K2 fiGould fc Horne 


helds for observin 


g. The 


2013 


Howell et ah 


2014 . 


WFIRST f Spergel et al.ll2015 h and Euclid I Penny et alJl2013 h Our method could also be 


applied in the process of held selection for LSST, if the G alactic disk 


frequently enough to allow detection of microlensing events (IGouldl 120131) . 


lelds are observed 
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Figure 7. Expected event rates for all bulge OGLE-III subfields with derived A^rc- The values are based on Eq.[3] 
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Name 

1 

[deg] 

b 

[deg] 
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4.708 
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